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^ ■ Abstract 



In Brans-Dicke theory the Universe becomes divided after inflation into many ex- 
ponentially large domains with different values of the effective gravitational constant. 
Such a process can be described by diffusion equations for the probability of finding 
a certain value of the inflaton and dilaton fields in a physical volume of the Universe. 
For a typical chaotic inflation potential, the solutions for the probability distribution 
never become stationary but grow forever towards larger values of the fields. We show 
here that a non-minimal conformal coupling of the inflaton to the curvature scalar, as 
well as radiative corrections to the effective potential, may provide a dynamical cutoff 
and generate stationary solutions. We also analyze the possibility of large nonper- 
turbative jumps of the fluctuating inflaton scalar field, which was recently revealed 
in the context of the Einstein theory. We find that in the Brans-Dicke theory the 
amplitude of such jumps is strongly suppressed. 
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1 Introduction 



After fifteen years of development of inflationary cosmology, the basic principles of this theory 
seem to be well understood, see e.g. Ref. [|TJ. However, there is still a very long way from these 
basic principles to the final theory. One of the main problems is the absence of the final version of 
the underlying particle theory. Fortunately, many properties of inflationary cosmology are very 
stable with respect to the change of its particular realization. In particular, most inflationary 
models predict a flat Universe (0=1) with a nearly scale-independent spectrum of density 
perturbations. Still some important deviations from the standard lore may appear when one goes 
from one theory of elementary particles to another. For example, in a certain class of theories, 
one can obtain an open or a closed || |j] Universe, or even a Universe consisting of different 
causally disconnected regions with all possible values of ft < 1 [Q. It would be important to 
find out whether some other modifications of the inflationary paradigm may appear when one 
implements it in different theories of fundamental interactions. 

An interesting playground for testing the robustness of various ideas about inflation is pro- 
vided by the Brans-Dicke theory ||. It remains to be seen whether this theory or some of its 
generalizations can have a sufficiently good motivation, e.g. from the point of view of string 
theory ||— 1|. In any case, however, some qualitatively new effects which appear in the BD 
inflationary cosmology may justify its investigation. 

One of these effects is the possibility to avoid the graceful exit problem of the first-order 
inflation |J. This possibility, however, requires the introduction of an effective potential and/or 
a nonminimal kinetic term for the Brans-Dicke field [|Hj], which makes the corresponding theory 
rather complicated and may lead to certain problems ||11|| . In some of the recent versions of these 
models the end of inflation occurs in the standard way, during the stage of slow rolling W2 . 



In what follows we will be interested in another specific effect which may appear in the 
inflationary Brans-Dicke cosmology. In this theory quantum fluctuations of the Brans-Dicke 
field during inflation driven by the inflaton field a lead to the division of the Universe into 
different exponentially large regions where the effective gravitational constant G = Mp 



2 



2nd 



can take all possible values from to oo, see Ref. [T3[ . 



This effect becomes especially interesting if one takes into account the process of self-repro- 
duction of inflationary universe. This process can be studied in a most adequate way by using the 
stochastic approach to inflation and investigating the probability distribution P p ((p,a,t), which 
gives the relative fraction of the volume of the Universe containing the inflaton field a and the 
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Brans-Dicke field at the moment t. A detailed study of the distribution P p was performed 
recently in Ref. [I3|, where it was shown that in many inflationary models based on the Einstein 



theory of gravity (i.e. the theory with a constant Brans-Dicke field) this probability distribution 
rapidly approaches a stationary regime. This means that if one takes a section of the Universe at a 
given time t and calculates the relative fraction of domains of the Universe with given properties, 
the result will not depend on the time t, both during inflation and after it. 

However, the stationary character of the probability distribution P p in general relativity is 
closely related to the existence of the Planck boundary, where the potential energy density V(a) 
becomes comparable with the Planck density Mp\ Typically the distribution P p (a,t) rapidly 
moves towards large a, for the reason that the volume of domains with large V(a) grows very 
fast. The distribution P p (a,t) becomes stabilized as it approaches the Planck boundary, where, 
as it is argued in [r4||-||16||, the process of self-reproduction of inflationary domains is less efficient, 



or inflation becomes impossible altogether. 



Meanwhile, in the Brans-Dicke theory the situation is more complicated ]T7|, |IBf . Since the 
effective Planck mass in this theory depends on <f), Mp(0) = — <ft 2 , the Planck boundary instead of 
being a point becomes a line 4 = -A V(a). Then the distribution P p (0, a, t) first moves towards 
the Planck boundary, and after that it slides along this boundary in the direction where the 
rate of expansion of the Universe becomes greater. Such runaway solutions typically lead to the 
distributions P p rapidly moving towards indefinitely large values of <fi and a. As a consequence, 
a typical value of the Planck mass at the end of inflation also becomes indefinitely large, and the 
effective gravitational constant vanishes. 



In fact, runaway solutions may appear even in the ordinary Einstein theory. This may happen 
if the inflaton field is nonminimally coupled to curvature due to the existence of the interaction 
— |£cr 2 -R, with £ < [16]. On the other hand, in the theories with £ > one may encounter many 



different regimes. Depending on the value of £, one may not have inflation at all, or one may have 
inflation without self-reproduction of inflationary domains, or inflation and self-reproduction with 
a stationary distribution P p . 

This suggests us to study a more general version of the Brans-Dicke theory, where not only 
the Brans-dicke field <p, but also the inflaton field a is nonminimally coupled to gravity. Some 



models of this type have been discussed before in [|T^, [2(|. Note that this framework seems to be 
much more natural than the standard one, where it is assumed that only one of the two scalar 
fields (the Brans-Dicke field) is nonminimally coupled to gravity. 

This paper is organized as follows. In Section [| we discuss the problem of runaway solutions 
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in the Brans-Dicke stochastic inflation. In Section ^ we derive a set of equations for inflation in 
the Brans-Dicke theory with a non minimally coupled inflaton scalar field. We study the case of a 
purely massive scalar field and find that it provides a natural cutoff for the motion of the inflaton 
and dilaton fields, resulting in the existence of stationary solutions to the diffusion of both fields. 
In order to test the stability of the results, we analyse the same model in the presence of a non zero 
self-coupling of the inflaton. We find that such a term is enough to destroy the stationarity of the 
probability distribution, making such a stationarity rather improbable. In Section [3J] we propose 
a different solution to the runaway behavior of the scalar fields by considering 1-loop corrections 
to the effective potential. For certain values of the parameters of the model, it is possible that 
the effective potential of the inflaton field acquires a maximum at exponentially large values of a. 
This provides a natural cutoff for the rate of inflation and makes the distribution P p stationary. 
Furthermore, the Planck mass in this model is exponentially larger than the only scale in the 
problem, the mass of the inflaton, thus naturally explaining a hierarchy of scales. In Section || 
we explore the possibility of very large quantum jumps of the scalar field [21] in the context of 
the Brans-Dicke theory. In Section [| we draw some conclusions. 



2 Runaway solutions in Brans-Dicke inflation 



V(a) 



(1) 



In this section we will introduce the problem of runaway solutions in Brans-Dicke cosmology. For 
a detailed analysis see Refs. |T^, . Consider the evolution of the inflaton field a with a generic 
chaotic potential in a JBD theory of gravity with dilaton field </>, 

S = [ cfixJ^q — ( 
J V y l8cu 

where Planck mass is written in terms of the dilaton field as Mp(0) = ^-(fi 2 - For generic inflaton 
potentials of the type V(a) = ^cr 2n , the equations of motion for the homogeneous fields in the 
slow-roll approximation are 

b n H 4> 2 
a 2 uj a 2 



H 



This implies that 



d 



+ -a' 

n 



(2) 
(3) 



which means that the fields move along a circle of radius r 2 = cf 2 + (p 2 , with (p 2 



Ref . |TB| . Inflation begins at the Planck boundary 



(j) ~ ^ a ^> a . This implies that if inflation begins at o~o 0o 
if we find a reason why o"o should be large, we will be able to explain why Mp 
is also large. 



2 2 

- cr , see 

Jf^b a n and ends somewhere at the line 

27r\/2n 

it will end at (f> ~ \/f c"o- Thus, 

2tt i 



C"0 
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This simple picture becomes much more complicated if one takes into account that in addition 
to the classical motion, the fields ef> and a experience quantum jumps with a typical amplitude 
during each time interval H~ l . One can easily verify that (for ^ < a ) these jumps are greater 
than the classical rolling of the fields for 

371-2 ^4 \ 4vr 2 4 

— 4 < V{a) < — 4 , (4) 
ur or 

the last inequality corresponding to the Planck boundary. Thus, if the field is not very far away 
from the Planck boundary (but still sufficiently far away so that our approach remains reliable), 
the motion of the field occurs not due to classical rolling but due to Brownian jumps in all 
possible directions. The jumps in the direction of greater Hubble constant lead to much more 
rapid expansion of space. As a result, very soon (in the synchronous time t) the main part of the 
physical volume of the Universe shifts towards the Planck boundary where the Hubble constant 
takes its greatest values. After that the field a can take all possible values along the Planck 
boundary, and therefore the value of the field <fi in the domains reaching the boundary of the end 
of inflation does not depend on the initial value a . 

However, the Hubble constant H(<p, a) is different in different parts of the Planck boundary. 
If it has a maximum at some point a = o~ max along the Planck boundary, then the main part of 
the volume of the Universe will be produced as a result of exponentially rapid inflation near this 
point. Therefore the main part of the volume of the Universe after inflation will be produced 
with M P ~ y/Jcw 0. 

Unfortunately, as we have shown in |17j this regime is not realized in the model ([I]) with 
V(a) = j^cr 2n , nor in the models with exponentially growing potentials. The reason is very 
simple: The Hubble parameter in such theories grows along the Planck boundary. Therefore 
the leading contribution to the volume of the Universe is given by the domains expanding with 
ever growing speed and containing indefinitely large values of the fields a and <fi. The diffusion 
equations describing this process have runaway solutions which describe distributions P p ((j),a,t) 
running towards infinitely large a. In many cases, the center of the distribution P p reaches 
infinitely large values of a within finite time. 

This does not necessarily mean that the corresponding theories are physically unacceptable. 
It may happen that our idea that we should live in a part of the Universe corresponding to a 
maximum of P p {4>, cr, t) at a given time t is incorrect; see Ref. [|l7j] for a discussion of this issue. In 



particular, as it was shown in JTJ], pp , the behavior of the distribution P p (4>, cr, t) depends on the 
choice of time parametrization. If one studies, e.g. the probability distribution P p (<f>, cr, r), where 
r = lna(t), its behavior may be quite different, and it may not exhibit any runaway solutions. 
Moreover, it is very easy to stabilize the distribution P p (4>, cr, t) by making the effective potential 
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V(a) very curved and unsuitable for inflation at large a . Still it would be interesting to study 
this phenomenon in a more detailed way and to see whether the distribution P p ((f),a,t) in the 
Brans-Dicke theory can be stabilized in a natural way. If this stabilization can occur only at very 
large values of the field a, then the typical value of the Brans-Dicke field <fi at the end of inflation 
will also be extremely large. This may give us a tentative explanation of the anomalously large 
value of the Planck mass. 



3 A more general theory 



Let us now turn to a scenario in which some of the problem discussed above could be resolved. 
For that purpose we will consider the classical evolution of the inflaton field with a generic chaotic 
potential, in the context of the Jordan-Brans-Dicke theory of gravity and a curvature coupled 
inflaton with a non minimal coupling £, 



S = J d Xy/—g 
In this theory the Planck mass takes the form 

16vr$ 



quj 2 2 2 



V{a) 



M 2 ( 



(5) 



(6) 



where $ plays the role of the Brans-Dicke scalar. The value of the Brans-Dicke parameter uo is 



bounded by the post-Newtonian experiments |]22j and primordial nucleosynthesis |23[ to be very 
large, to > 500, and therefore it is reasonable to use the approximation uj ^> 1 in the following 
analysis. The parameter £, on the other hand, is unconstrained and could take positive or negative 
values. Here we will consider the case of a small and negative £. 



The equations of motion for the theory (|5|) can be written as 

W = -V'(a)-£aR , 

We can then write the exact equations for the homogeneous fields in flat space (k = 0) as 



+ v 2 2 + (l-6OVV 
2oo ' 



AV(a) - aV'(a) 



(7) 



(8) 



(9) 
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—(l — 



2u 



i? = W(a) - QiaV\a) - (1 - 6£)<x 2 - (l + — 



i-0 2 -^ 2 ) 3tf 2 = K(a) + 6^ ( xa + ia 2 -|-^+i0 2 



(10) 



where R = 12H 2 + 6H. Note that for £ = and $ = <J) 2 /8uj we recover the usual BD equations. 
From now on we will drop the last two terms of Eqs. ( |i~0"D since they are subleading for large u;. 



During inflation, we can write the equations of motion of the homogeneous fields <fi and a for 
small |£| and large a;, in the slow-roll approximation (0 <C iJ0 <C H 2 <p), as 

• #0 



<7 



UJ 



1 



41/ a 



+ --5 



(7- 



11' 



# 2 



4\/(cr) 
3 </> 2 + /x cr 2 



where /x = 4a>|£|. Note that for small \x they reduce to the usual BD equations. 



3.1 Massive curvature-coupled inflaton 

Let us consider in some detail the theory fl5|) with V(a) = ^m 2 a 2 at large a and very small £ < 0. 
In this case it is useful to make a change of variables to polar coordinates r 2 = <ft 2 + (p 2 , z = 
<£> 2 = 2(l + 6|£|)cr 2 ~ 2a 2 . In these coordinates, the equations of motion (|Tl|) take the form 

• ZH (z 2 + 1-^/2) , 



rH /j/2 



uj 1 + z 2 ' 
if 2 = _ m2 



3 z 2 + /i/2 

One can solve these equations parametrically in the z-variable, for /i ^ 2, 
c?r r/i/2 



z(l + z 2 )(z 2 + l-///2) ' r ° r( ^ o) 

r(z) = r ' 1 



which gives 



0(2) = r(z) (1 + z 
ip(z) = r(z) (1 + z 



-2N-1/2 



2\-l/2 



(12) 



(13) 



(14) 
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For |£| <C 1, the condition for existence of inflation (\H\ < H 2 ) is a 2 < V, 

(z 2 -^/2) 2 <^(z 2 + f ,/2), (15) 

which gives a slightly modified condition for the end of inflation, z 2 + ft/2 = 3u/2. 

For very small /i, the trajectory is very close to a circle in the plane (</?,</>). However, for 
arbitrary /i the scalar fields move along a more complicated trajectory. There are regions of 
parameter space for which inflation never ends. For instance, if // > 2 and z 2 < (/i — 2)/2, then 
i < and inflation never ends (unless the field a tunnels (diffuses) to its large values quantum 
mechanically). This is a rather surprising result since there is nothing singular in our equations 
(|12|) at fi = 2. The reason of this unexpected behavior is that the field a in the regime -C a for 
fj, > 2 grows much faster than the field 0, so it never approaches the end of inflation with 3> a. 

The value of the effective Planck mass at the end of inflation is given by (see Eq.[T3|) , 

/2vr . . fz 2 + u/2\ 1/2 [2n ( 2 - 

where we have used the fact that z e ^> 1 and z a <C 1. Note that the value of the Planck mass at 
the end of inflation is very sensitive to initial conditions and, for the reason explained above, it 
becomes exponentially large when fi — >• 2: Mp = ^Atx\^\ (fi exp (z~ 2 ). 

We should now discuss possible initial conditions and the regime of self-reproduction of the 
Universe in this model. The Planck boundary in this theory looks like an ellipse 

\ 2 2 
l2 i , , / „ °max \ A^max 

+ fj, [a 



2 

C"m ax — — 



(17) 



•57T 



There is no constraint on fi; however, we will assume that £ and \ju) are both of the same order 
of magnitude. We will see that this is a consistent approximation in this scenario. 

The amplitude of quantum fluctuations of the scalar fields, 5(f) and 5a, whose wavelengths 
are stretched beyond the horizon and act on the quasi-homogeneous background fields like a 



stochastic force, can be computed as in |17], [L8j and they are given by (for u ^> 1 and |£| <C 1) 



H r H 

, OCT ~ . 

2vr ' 2vr 



Density perturbations in our model can be calculated as in Refs. [ 17] , |18| , 2^| . For the theory 



8 



V(o~) = ^m 2 a 2 , in the large uj limitf], we find 

Sp 50 m 



(19) 



p M P 

Note that the larger is the Planck mass at the end of inflation in a given region of the Universe, 
the smaller will be density perturbations in this region. This suggests that the large value of 
Planck mass Mp in our part of the Universe may be related to the small value of the amplitude 



p 

Sp 

M P L p 



of density perturbations -j ~ 5 • 10 JT^j- In other words, instead of two independent small 



parameters, 4j- <C 1 and — ~ 5 ■ 10 5 , we have only one. 



As we already discussed in the previous section, those inflationary domains where quantum 
jumps during the time if -1 are more important than their classical motion, enter the regime of 
self reproduction. One can easily show that in our case the region where this regime is possible 
is the interior of the ellipse 



(2 



+ V[(T- — 



T \ 2 Vol 



21 4 (20) 



In our model the Hubble constant takes its maximal values along the line = (i.e. z — 0). 
Along this line there is a wide region where self-reproduction is possible. This region begins at 
the Planck boundary flT?D and ends at the self-reproduction boundary (p0|): 

mm , . 

<^< : y^,,,, • (21) 



8tiV2|£| 4vrv / 6|e| 3/2 

Since this interval is limited, runaway solutions here are impossible. In such a situation the 
probability distribution P p (<f), o~, t) should be stationary, with a maximum somewhere at this line, 
in the interval (pl|). 

Note, however, that this maximum should be relatively smooth, since the value of the Hubble 
constant does not depend on a along the line = 0: 

m 

H(d> = 0) = -= . (22) 

For a complete investigation of the probability distribution P p ((p,a,t) and of the value of Mp at 
the end of inflation one should solve numerically the diffusion equations for P p (<f), a, t), as we did 



in Ref. [17]. To get a rough idea of the resulting distribution of possible values of Mp at the end 



3 Note that during the last stages of inflation, there is an approximate equivalence of the Einstein and Jordan 
frames. 
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of inflation one can simply take all points at the boundary of self-reproduction and treat them as 
initial conditions for the solutions (|IBD, (P^[). However, one can get a much better picture if one 
takes into account that the main contribution to the probability distribution on this boundary 
is given by the part of the boundary of the region of self-reproduction with the highest value of 
H. This value is achieved for = 0* (|20|) at the line = 0. The peak cannot be wider than 



^(J* in the 0-direction, and in fact it is expected to be much more narrow. Estimates of the 
values of the effective Planck mass generated after inflation suggest that Mp can take all possible 
values from to 00, but typical values are much greater than -j=g. If, for example, one considers 
the evolution of domains with initial values of the fields a ~ a* and > ^- (due to unavoidable 



quantum jumps near = 0), Eq. (0) yields 



> J3^L . (23) 



M P ~ v r \2- f i / 

One can then estimate the allowed values of |£| from the amplitude of density perturbations 
(|T5|) . For n ~ 1 we find |£| < 4 ■ 10~ 4 , which is consistent with the constraint uj > 500 and our 
condition fx ~ 1. Note also that ~ \/&n exp ( — 2j|p) f° r Z 1 — > 2, and it becomes very easy to 
obtain a very small ratio of ^-p- even for not too small values of |£|. Therefore, we have a natural 
realization of our model without the need for very small numbers. In fact one could argue that 
uj ~ 500 is not such a large number, if we understand it as the dimensionless coupling (3 of matter 
to the dilaton in the Einstein frame p| , 2f3 = {2uo + 3)~ 1/2 ~ 0.03. 



3.2 Self-coupled inflaton 

We must now consider possible corrections to the inflaton potential, in order to see whether 
our results are stable with respect to small modifications of the theory. Let us add an inflaton 
self-coupling to the theory, V(a) = ^m 2 a 2 + |cr 4 . The corresponding equations of motion read 

. zH 2m 2 (l + z 2 )(2z 2 + 2 - //) + Ar 2 (l + 2^ 2 ) 



to 4m 2 (1 + z 2 ) + Ar 2 

. _rH 2m 2 (l + z 2 )fi - Xr 2 z 2 

r = ~ (4m 2 (l + z 2 ) + Ar 2 )(l + z 2 ) ' ^ ' 

9 uj 4m 2 (l + z 2 ) + Ar 2 
H z - 



6 (1 + z 2 ){2z 2 + fj) 



There is a bifurcation line in the plane (z, r) that separates domains for which inflation never 
ends. This is the line z = 0, 

,2_^~2 Ar 2 l + 2z 2 

2 4m 2 1 + ^ 2 ' 1 ' 
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For fi > 2 and r 2 < 2m 2 (/i — 2) /A, there are values of z for which z < 0, which could suggest 

l 

2 



just as in the previous model with V(a) = \w?a 2 , that the end of inflation cannot be reached. 



However, as r increases ([24]), z eventually changes sign and inflation ends. 
The boundary of self-reproduction in this theory is given by 

2 + ua A ( 1 _f\ = vajtG (l + ^l) , (26) 



a 



where v and ex* were defined in (f2"UD. 

The Planck boundary for this theory is 

2 + /icr 2 = /icr max cr (1 + ^—1 . (27) 




where <7 max was given in ( |T?j ) and cr c = 2m 2 / v^A- At small a this boundary looks like an ellipse, 
but at large a it becomes a straight line. For cr max < o c these two parts of the Planck boundary 
are disconnected: both lines cross the a axis at different points. For cr max > a c these two lines 
form a continuous curve. The greatest values of the Hubble parameter appear along the part 
of the Planck boundary with a —>■ oo. Therefore in this model there always exist runaway 
solutions describing the probability distribution P p (<j), cr, t) rapidly moving towards indefinitely 
large values of a and 0. In other words, the existence of a stationary regime which we have found 
in the previous model appears to be unstable with respect to a small modification of the effective 
potential of the inflaton field. It makes this model less attractive. 



4 One loop corrections 

Now let us return back to the simple Brans-Dicke theory (|l|) with the inflaton field minimally 
coupled to gravity, but let us take into account quantum corrections to the effective potential 
V(a). In the one- loop approximation one can represent the effective potential in the following 

way: 

V(a) = j(a 2 -a 2 ) 2 + (3aHn- , (28) 

where a\ = m 2 /X is the value of the inflaton field at the minimum of the potential, and (3 depends 
on the values of coupling constants in the theory. For example, if the scalar field a interacts with 
the fermion field with the coupling constant h, one obtains j3 = (9A 2 — Ah 2 ) /32tt 2 ||26|| . Note, 
that (3 will be negative for h 2 ^> A 2 . In that case, a new maximum will appear at exponentially 
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large a, 

Omax ^ (T exp — — > <7 , (29) 

followed by a very sharp fall-off to negative values of the potential, rendering the vacuum unsta- 
ble ||26|| . In our case this sharp fall-off may act effectively as a natural boundary for the diffusion of 
the inflaton field a. In such a situation the probability distribution P p (<p, a, t) becomes stationary, 
with the maximum concentrated near the Planck boundary at o = o" max - 

In the theory ([28|), the equations of motion for a ^> o~o an( f £ = 0, become 

• H<f> 



V'(a) H4> 2 



H 



3H ujo 
2 AooV(a) 



(30) 



For do <C o <C o" max the classical motion is circular, just as in the theory Act 4 , and we can define 
polar coordinates (r,z) as usual, with z = <p/u. Therefore one can obtain the following estimate 
of the Planck mass after the end of inflation: 

M P ~y-a max ^-a exp— . (31) 

Note that in this case the effective Planck mass for \/3\ A becomes exponentially large. In the 
regime when it is much greater than <7o the amplitude of density perturbations is given by the 
standard expression ^ ~ 10 2 \/A. Thus, one does not get anything new from the point of view 
of density perturbations, but one can obtain a natural explanation of the very large value of the 
Planck mass as compared with other mass scales in the theory. 



5 Non-perturbative effects 

All results obtained in the present paper are related to the probability distribution P p ((p,a,t), 
which shows the fraction of the volume of the Universe with the fields <fi and a at a given moment 
of time t in synchronous coordinates. We should emphasize again that this is not a unique choice 
of measure in quantum cosmology. One could study, for example, the probability distribution 
P p (0, a, r), where r ~ lna(t). This distribution is considerably different from P p ((f),a,t). This 
does not allow us to make unambiguous predictions until the issue of measure in quantum cosmol- 
ogy is resolved, see a discussion of this issue in Refs. ]14[ [17[ [L6[ . Still we believe that investigation 



of P p ((j),a,t) gives us very interesting information about the structure of inflationary universe. 
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Recently it was shown, in the context of chaotic inflation based on the usual Einstein theory, 
that the main fraction of volume of the Universe in a state with a given density p at any given 
moment of time t is concentrated near the centers of deep exponentially wide spherically sym- 
metric holes in the density distribution PT|. F° r t ne reason discussed above, interpretation of 



this result is not unambiguous, and we are not sure that it implies that we must live near the 
center of a spherically symmetric void. However, we think that this result is very nontrivial and 



deserves further investigation [27 



Here we would like to study whether a similar effect occurs in the context of the Brans-Dicke 
theory. As we will see, the effect does take place, but its amplitude is significantly different. In 
order to understand it, we should briefly remind the origin of this effect in the Einstein theory, 
with one scalar field a (the inflation field), and then we will make a generalization to the Brans- 
Dicke case. 

The best way to examine this scenario is to investigate the probability distribution P p (a,t). 
The distribution P p (a,t) obeys the following diffusion equation (see Ref. jnj and references 
therein) : 

dP p 1 d fH 3 / 2 (a) d /# 3/2 (a) d n V'(a) 



Here we temporarily use the system of units Mp = 1. One may try to obtain solutions of equation 
(|32|) in the form of the series P P (<J, t) = J2^L\ eXst 7Ts(c)- In the limit of large time t only the term 
with the largest eigenvalue Ai survives, P p (a,t) = e Xlt ni(a). The function 7Ti in the limit t — > oo 
has a meaning of a normalized time-independent probability distribution to find a given field a 
in a unit physical volume, whereas the function e Alt shows the overall growth of the volume of 
all parts of the Universe, which does not depend on a in the limit t — ► oo. In this limit one can 
write Eq. fl32|) in the form 



1 d f H^(a) d ( HW(a) \\ d ( V'(a) , ,\ , x ( , ^ 



In the simplest theory with V(a) = fa 4 and H = y^a 2 , Eq. © reads |4[: 

// // 6 9\ / 6 15 367r Ai /67T\3/2\ 
1 ^Aa 5 o> Ucr 6 a 2 Act 4 trt 6 V A ' ' 

This equation can be solved both analytically and numerically. The result is that the eigenvalue Ai 
is given by d(X) H max . Here d(X) is the fractal dimension, which approaches 3 in the limit A — > 0, 
while H max is the maximum possible value of the Hubble constant during inflation, which in our 
case corresponds to its value at the Planck boundary V(a) = Mp = 1. This gives H max = 2\f^f- 
Thus, in the small A limit one has Ai = 3if max = 2\^6n w 8.68. 
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The distribution tt\ depends on a very sharply. One can easily check that at small a the 
leading terms in Eq. (^) are the second and the last ones. (This means, in particular, that the 
diffusion terms in Eqs. fl32|) , (|33l) can be neglected.) Therefore the solution of equation ([34] ) for 
small a corresponding to the last stages of inflation is given by 



7Ti ~ ffVT A i . (35) 
This is an extremely strong dependence. For example, ix\ ~ a 10 * for the realistic value A ~ 10~ 13 . 

Consider all inflationary domains which contain a given field a at a given moment of time t. 
Let us try to find a typical value of this field in those domains at the earlier moment t — H~ l . 
In order to do it, one should add to a the value of its classical drift aH~ x . One should also add 
the amplitude of quantum jumps Act. The usual estimate of the magnitude of a typical jump 
is ±^-. This is a correct estimate if we are interested in a typical amplitude of jumps at any 
given point. However, if we are considering all domains with a given a and trying to find all 
those domains from which the field a could originate, the answer for may be quite different. The 
total volume of all domains with a given field o at any moment of time t strongly depends on 
cr: Pp(cr) ~ tti(o") ~ crV / ^" Al , see Eq. (|35|). This means that the total volume of all domains 
which could jump towards the given field a from the value a + Act will be enhanced by a large 

additional factor ~ [ \ + ^-j x \ On the other hand, the probability of large jumps 

Act is suppressed by the Gaussian factor exp(— 27r ^ g2 )- The product of these two factors has a 
sharp maximum at Acr = Aicr • In other words, most of the domains of a given field a are 
formed due to the jumps which have definite sign (they are decreasing the value of the scalar 
field), and which are greater than the "typical" ones by the amplification factor iV = Aid. In the 
usual notation, this amplification factor is given by 

" = < 36 » 

where Ai = 3if max = 2\^6nMp « 8.68Mp. For a ~ 4.5Mp (the scale at which the large scale 
structure of our part of the Universe has been formed), the amplification factor in our theory is 
about 40. However, the value of this factor is very sensitive to our assumptions concerning the 



Planck boundary, it can be much bigger or much smaller than 40. As explained in ||21|| , this effect 
does not alter the standard theory of density perturbations in inflationary universe, but it puts 
these perturbations on the top (or, more precisely, to the bottom) of the distribution of the scalar 
field a which appears as a result of its large jumps. 

We can now return to the Brans-Dicke theory with u ^> 1. In this case both fields and 
a move and fluctuate, and therefore one should write a two-dimensional diffusion equation for 
these fields W7\ . However, at the last stages of inflation the classical motion of the scalar field <j) is 
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very slow, and, as we argued above, diffusion in the first approximation can be neglected. Then 
the problem reduces to the one we have already solved, and the amplification coefficient will be 
given by eq. fl36l). The only difference is that now the value of the coefficient Ai = 3if max will be 
determined not by the Planck mass at the end of inflation, but by the much smaller Planck mass 
at the place near the Planck boundary corresponding to the peak of the distribution P p (cj),a). 
If the distribution is stationary due to the existence of some kind of boundary at a = cr max , 
then the peak is concentra ted at the Planck boundary near a max . In this case one can show that 
Ai = 3H(a max ) = 2<7 max \/ 37rV% for V(cr) = | cr 4 . Meanwhile, the typical Planck mass after 
inflation is given by Mp ~ y^Omax ||17|| . This leads to the following realization of Eq. for 
the Brans-Dicke theory: 



which gives N ~ 4.5^^, for cr ~ 4.5Mp. Note that for A ~ 10 13 and uj ~ 10 3 the factor 
N becomes very small, and therefore all the nonperturbative effects discussed above will be 
negligible, contrary to the case in general relativity |?T] . 



Note that in this investigation we have assumed the existence of the upper boundary cr max . 
The situation will be quite different if there were no stationary solutions for P p . In such a case 
the nonperturbative effects will be extremely strong. However, then we will meet the problem of 
runaway solutions which would suggest that Mp — > oo in the main part of the physical volume 
of the Universe. In this paper we have shown that in principle it is possible to avoid runaway 
solutions in the Brans-Dicke theory. But then from our results it follows that the same trick 
which makes the distribution stationary simultaneously kills the nonperturbative effects in the 
inflationary Brans-Dicke cosmology. 



6 Conclusions 



Investigation of the probability distribution P p gives us a lot of interesting information about the 
properties of the inflationary universe. Some of these properties (such as the very existence of the 
regime of self-reproduction and the fractal structure of the Universe) do not depend on the choice 
of time parametrization, and therefore their interpretation is relatively straightforward. Some 
other properties of P p do depend on the choice of time parametrization. Sometimes it is not 
enough to know P p , we need to know also whether all parts of the Universe which are described 
by this distribution are equally well suited for existence of life. In such situations interpretation 
of the results becomes increasingly speculative. Nevertheless we believe that even in these cases 
investigation of P p and attempts of its interpretation can be very useful. We are learning how to 
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formulate questions in the context of quantum cosmology. In some cases we are obtaining results 
which look obviously incorrect or contradict observational data; then we may conclude that we 
are using quantum cosmology in a wrong way. In some other cases quantum cosmology allows us 
to obtain important results which cannot be obtained by other methods. This may be considered 
as an indication that we are on the right track. Thus, by this trial and error method we may 
finally learn how to use quantum cosmology. 

To give a particular example one may consider the old issue of the wave function of the 
Universe. The two most popular candidates are the Hartle-Hawking wave function, which in 

i6V W J t@' anc ^ ^ e tunneling wave function 
— igyfe j p9| . Neither of these two functions was rigorously derived. For some reason which 
is not related to its derivation, the square of the Hartle-Hawking wave function correctly describes 
tunneling between two different de Sitter universes with different values of V(<f)) [[30 1. However 



if one makes an attempt to apply it to the probability of creation of the inflationary universe, 
one comes to a physically incorrect conclusion that it is much easier to create an infinitely large 
inflationary universe with V(4>) — > than a Planck-size universe with V(4>) ~ Mp. The tunneling 
wave function leads to a qualitatively correct description of quantum creation of the Universe 



"from nothing," but one should not uncritically apply it to, e.g. the formation of black holes [31 



The reasons for the limited applicability of each of these functions are explained in Ref. [|TJ 
whereas in Ref. |14| it is shown that expressions of the type of exp ( Jyfej ) or exp(— 
appear in many problems of quantum cosmology which are not related in any obvious way to the 
original "derivations" of the Hartle-Hawking and tunneling wave functions. 

Something similar may occur with our investigation of the probability distribution P p . This 
distribution is certainly very useful, but sometimes it becomes tempting to use it in the situations 
where its interpretation is ambiguous and the final success is not guaranteed. However, the 
possibility to obtain very strong results and to look at the interplay between particle physics 
and cosmology from an entirely new point of view suggests us to continue this trial-and-error 
investigation. Here we may mention our attempt to address the cosmological constant problem 



in the context of the Starobinsky model []T6[ and the possibility to explain why in the main part 
of the Universe the scalar-tensor theories of gravity are reduced to the Einstein theory In 
this paper (see also our previous publications []T3|, [T7|, [TB[) we demonstrated that in the context 



of inflationary Brans-Dicke theory it may be the possible to explain the anomalously large value 
of the Planck mass Mp. We have shown also that the structure of our part of the Universe which 



could appear as a result of nonperturbative effects in quantum cosmology [51] may be extremely 
sensitive to the properties of the theory at nearly Planckian densities and to the presence or 
absence of stationary solutions for the distribution P p . At the very least, what we have found 
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can be considered as a description of rather nontrivial properties of hypersurfaces of a given 
synchronous time in the inflationary universe. However, we hope that some of our results may 
have deeper physical significance. 
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